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ABSTRACT 

We present spectroscopic and photometric observations of the Type Iln supernova 
(SN) 2008iy. SN 2008iy showed an unprecedentedly long rise time of ~400 days, making- 
it the first known SN to take significantly longer than 100 days to reach peak optical 
luminosity. The peak absolute magnitude of SN 2008iy was M r w —19.1 mag, and the 
total radiated energy over the first ~700 days was ^2 xlO 50 erg. Spectroscopically, 
SN 2008iy is very similar to the Type Iln SN 1988Z at late times, and, like SN 1988Z, 
it is a luminous X-ray source (both supcrnovae had an X-ray luminosity Lx > 10 41 
erg s _1 ). SN 2008iy has a growing near-infrared excess at late times similar to several 
other SNe Iln. The Ha emission-line profile of SN 2008iy shows a narrow P Cygni 
absorption component, implying a pre-SN wind speed of ^100 km s _1 . We argue that 
the luminosity of SN 2008iy is powered via the interaction of the SN ejecta with a 
dense, clumpy circumstellar medium. The ~400 day rise time can be understood if the 
number density of clumps increases with distance over a radius ~1.7 xlO 16 cm from 
the progenitor. This scenario is possible if the progenitor experienced an episodic phase 
of enhanced mass loss < 1 century prior to explosion or if the progenitor wind speed 
increased during the decades before core collapse. We favour the former scenario, which 
is reminiscent of the eruptive mass-loss episodes observed for luminous blue variable 
(LBV) stars. The progenitor wind speed and increased mass-loss rates serve as further 
evidence that at least some, and perhaps all, Type Iln supernovae experience LBV-likc 
eruptions shortly before core collapse. We also discuss the host galaxy of SN 2008iy, 
a subluminous dwarf galaxy, and offer a few reasons why the recent suggestion that 
unusual, luminous supernovae preferentially occur in dwarf galaxies may be the result 
of observational biases. 
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1 INTRODUCTION 

The recent development of synoptic, wide-field imaging has 
revealed an unexpected diversity of transient phenomena. 
One such example is the discovery of a new subclass of very 



luminous supernovae (VLSNe; e.g., Ofek et al. |2007 Smith 


et al. 2007 Quimby et al. 


2007b 


. While events of this 


nature are rare (Miller et al. 


2009a 


Quimby et al. |2009l, 



each new discovery serves to bracket our understanding of 
the physical origin of well-established classes of SNe and 
does, in principle, demand an increased clarity in our un- 
derstanding of the post-main-sequence evolution of massive 
stars. 



* E-mail: amiller@astro.berkeley.edu. 
f Sloan Research Fellow. 



Yet another unusual transient was discovered by 
the Catalina Real-Time Transient Survey (CRTS; |Drake| 
et al. 2009a I, which they announced via an ATel as 
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CSS080928: 160837+041626 on 2008 Oct. 07 UlQ ( [Drake et| Table x . Deep sky observations of SN 2008iy. 



al. 



2008). The transient was classified as a Type Iln SN 
2009 Mar 



with a spectrum taken on 2009 Mar. 27 (Mahabal et al. 
[2009{ see |Schlegel|[l996l for a definition of the SN Iln sub 
class, and Filippcnko 1997 for a review of the spectral prop- 
erties of SNe). The SN was later given the IAU designation 
SN 2008iy (jCatelan et ~aI7~| [2009]) . |Mahabal et~al~| p009| 



noted that the transient was present on CRTS images dat- 
ing back to 2007 Sep. 13; however, it went undetected by 
the CRTS automated transient detection software until 2008 
because it was blended with a non-saturated, nearby (~ll" 
separation) star. 

Here, we present our observations and analysis of 



SN 20 08iy, which peaked around 2008 Oct. 29 ( |Catelan et al 
| |2009 1 and had a rise time of ~400 days. This implies that 



SN 2008iy took longer to reach peak optical brightness than 
any other known SN. Type II SN rise times are typically < 
1 week (e.g., S Ne 2004 et and 2006bp ,|Li et al.|2005||Quimby 
|et al. |[2057a] see also |Patat et al~||1993| |Li et al. ||2010a| 7 
and have never previously been observed to rise > 100 days, 
let alone 400, making SN 2008iy another rare example of 
the possible outcomes for the end of the stellar life cycle. 
In addition to an extreme rise time, SN 2008iy is of great 
interest because the unique circumstellar medium (CSM) in 
which it exploded may provide a link to very long-lived SNe, 
such as SN 1988Z, and thus provide clues into the nature of 
their progenitors. 

This paper is organised as follows. Section 2 presents 
the observations. The data are analysed in Section 3, and the 
results are discussed in Section 4. We give our conclusions, 
as well as predictions for the future behaviour of SN 2008iy, 
in Section 5. 



2 OBSERVATIONS 
2.1 Photometry 

The field of SN 2008iy, which is located at a = 
16' 1 08 m 37.27 s , S = +04°16'26.5" (J2000), was imaged multi- 
ple times by the Palomar Quest survey, and those data have 
been reprocessed as part of the Deep Sky project]^] (DS 



Nu- 



gent|20 09). The best constraints on the explosion date come 
from DS imaging: in a coadd of two images from 2007 Jul. 
05, we do not detect the SN down to a 3cr limit of i > 20.9 
mag. DS images are best approximated by the Sloan Digital 
Sky Survey (SPSS; [Adeiman- McCarthy et al. 12008] ) i-band 
filter. The field of SN 2008iy has well-calibrated SDSS pho- 
tometry which we use to calibrate the DS images. The ob- 
served DS magnitudes are shown in Fig. [IJand summarised 
in Tabled] 

To constrain the rise time of SN 2008iy we visually esti- 
mate the possible range of magnitudes for the SN from im- 
ages on the CRTS websit^] based on a comparison to SDSS 
images. As previously mentioned, SN 2008iy is blended with 
a nearby star and so the automated photometry produced 
by the CRTS does not detect the SN on these epochs. On 



date 


mad a 




(MJD) 






54618.35 


17.37 


0.12 


54621.35 


17.38 


0.12 


54623.33 


17.35 


0.12 


54626.32 


17.34 


0.12 


54627.36 


17.37 


0.12 


54628.29 


17.37 


0.12 


54629.30 


17.36 


0.12 


54638.21 


17.34 


0.12 


54651.20 


17.28 


0.12 



a Observed value; not corrected for Galactic extinction. 

^ The calibration uncertainty dominates over the statistical 

uncertainty with values of ~0.12 and ^0.01 mag, respectively. 

2007 Sep. 13 the SN was between 18.3 and 19.0 mag, while 
on 2008 Feb. 19 the SN was observed between 17.9 and 18.3 
mag. 

Near-infrared (NIR) observations of SN 2008iy were 
conducted with the 1.3-m Peters Automated Infrared Imag- 
ing Telescope (PAIRITEL; |Bloom et aT~||2006[ ) starting on 
2009 Apr. 13. PAIRITEL observes simultaneously in the J, 
H, and K s bands. Observations were scheduled and executed 
via a robotic system, and the data were reduced by an auto- 
mated pipeline ( jBloom et al. |2006[ ). The SN flux was mea- 



sured via aperture photometry using SExtractor (Bertin & 
Arnouts||1996 1, calibrated against the Two Micron All Sky 
Survey (2MASS [Skrutskie et al. |2006[ ). Filtered photometry 
of SN 2008iy is shown in Fig. [2] and summarised in Table [2] 
Ground-based optical observations of SN 2008iy were 
obtained using the 1.0-m Nickel telescope located at Lick 
Observatory (Mt. Hamilton, CA, USA) starting on 2009 
Apr. 18. BVRI photometry f rom the Nicke l was measured 
using the DAOphot package ( |Stetson|[l987[ ) in IRAFQ and 
transformed into the Johnson-Cousins system. Calibrations 
for the field were obtained on ten photometric nights with 
the Nickel telescope. The Nickel photometry is shown in 
Fig. [2] and summarised in Table [3] 

SN 2008iy was observed by the space-based Swift ob- 
servatory on 2009 Apr. 15 and 16. Swift observed the SN 
simultaneously with both the X-ray Telescope (XRT 



rows et al. 



Bur- 



2005 I and the Ultraviolet / Optical Telescope 
(UVOT; |Roming et al.~||2005| ). We downloaded the Level-2 
UVOT images, and measured the U, B, and V photometry 



using the recipe of |Li et al. | ( |2006[ ). The UV filters (UVM1, 
UVW1, and UVW2) were calibrated using the method of 
Poole et al~| p008). Our final UVOT photometry is sum- 
marised in Table |4| The XRT observed the field for a to- 
tal of 10.71 ks. We extract the 0.3-10.0 keV counts from 
an extraction region of 64 pixels (~ 2.5'), where we fit the 



point-spread function model (see Butler fc Kocevski||2007| 
at the centroid of the X-ray emission. From the native XRT 



1 UT dates are used throughout this paper unless otherwise 
noted. 

2 http:/ /supernova. lbl.gov/~nugent/deepsky.html. 

3 http://crts.caltech.edu/. 



4 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation. 
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Figure 1. Left: Apparent optical light curve of SN 2008iy showing the long rise time of the SN, including data from Deep Sky (this 
work), the Nickel telescope (this work) , [Catelan et al, | ( |2009| l, and the CRTS website. Note: conservatively, we take the start time of 
the SN (t = 0) to coincide with the first epoch where CRTS detects the source; however, the true explosion date is likely prior to 2007 
Sep. 13, between the last Deep Sky non-detection and the first CRTS detection. Notice that the light curve is very broad, and almost 
symmetric about the peak. For comparison we show the light curve of the standard Type II-P SN 1999em (data from |Leonard et al.| 
1|2002[l and the Typ e Iln SN 1988Z, which shows many similarities to SN 2008iy at late times (see Section 3.2; SN 1988Z data from 
Turatto et al. |1993[ l, as they would have appeared at the redshift of SN 2008iy. Right: Three panels with images of the field of SN 2008iy. 
Each image is ~ 2.5' X 4', with north up and east to the left. Black crosshairs mark the location of the SN. The bottom image, marked 
(c), shows the last non-detection from the Deep Sky data. The middle (b) and top (a) images show the CRTS detections on day and 
day 153, respectively; SN 2008iy can clearly be seen in each of them. The SN was not flagged as a transient in either the middle or top 
image by the CRTS software. 



astrometry the centroid of the X-ray emission, located at a 
= 16 h 08 m 37.23 s , S = +04°16'26.7" (J2000), with a radial 
uncertainty of 5" (90% confidence interval), coincides with 
the optical position of SN 2008iy, suggesting that the X-rays 
are from the SN. In an astrometric fit relative to 2MASS we 
measure the position of the star ~11" from the SN to be 
a = 16 h 08 m 36.98 s , S = +04°16'16.4" (J2000), with a root- 
mean square scatter of 0.19" in a and 0.18" in 5 for our 
astrometric solution. This star is therefore well outside the 
90% confidence interval for the location of the X-ray emis- 
sion. Furthermore, there is no catalogued X-ray source at 
the position of the SN in the ROSAT catalogue ( |Voges et] 
al. ||1999 |, also suggesting that the star located ~11 from 
the SN is not the source of X-ray emission. We measure a 
background-subtracted 0.3—10.0 keV count rate of (19 ± 6) 
x 10~ 4 counts s -1 . 



side of the double-arm (blue+red) LRIS system. The spec- 
tra were reduced and calibrated using standard procedures 
(e.g., Matheson et al. 20001. Clouds were present on the 



night of 2009 Apr. 25, and the observing setup from that 
night did not include any wavelength overlap between the 
red and blue arms of the Kast spectrograph; thus, the ab- 
solute flux calibration on this night is less certain than on 
the other nights. For the Kast spectra, we observed the SN 
with the slit placed at the parallactic angle, so the relative 
spectral shapes should be accurate ( Filippenko|19 82 ). LRIS 
is equipped with an atmospheric dispersion corrector; nev- 
ertheless, we observed the SN at the parallactic angle on 
this night as well. In Fig. [3] we show the spectral sequence 
of SN 2008iy, and a summary of our observations is given in 
Tabled 



2.2 Spectroscopy 

Spectra of SN 2008iy were obtained on 2009 Apr. 18, 2009 
Apr. 25, and 2009 Jul. 23 with the Kast spectrograph on 
the Lick 3-m Shane telescope ( Miller fc Stone|1993 |. An ad- 
ditional spectrum was obtained on 2009 Sep. 22 with the 



Low-Resolution Imaging Spectrometer (LRIS; Oke et al. 
1995 1 on the 10-m Keck I telescope. To improve the spectral 



resolution at the location of the Ha emission feature from 
SN 2008iy, we used the 1200 lines mm -1 grating on the red 



2.3 Host Galaxy 

There is no catalogued redshift for the host galaxy of 
SN 2008iy. To determine the redshift of the SN we exam- 
ined the LRIS spectrum for the presence of narrow emission 
lines. Based on a fit to the [O III] AA4959, 5007 doublet we 
adopt a heliocentric redshift of z = 0.0411 for SN 2008iy (not 
corrected for peculiar motions or Virgo infall). This corre- 
sponds to a luminosity distance of dh = 179 Mpc (Ho = 71 
km s _1 Mpc _1 , Q,m = 0.27, Qa = 0.73), which we adopt 
throughout the remainder of this paper. 
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Figure 3. Spectral sequence of SN 2008iy (shown in black) and a single, late-time spectrum of SN 1988Z from our spectral database 
(shown in red), which was observed 504 days after discovery. To the left of the spectra we give the time since discovery, in rest-frame 
days, at which each spectrum was taken. Prominent spectral features of SN 2008iy are labelled at the top of the figure. The spectra show 
little evolution between days 560 and 711. Fe II multiplets 42, 48, and 49 can also be seen, though we have not labelled those lines. The 
spectra have been corrected for Galactic reddening (E(B — V) = 0.065 mag; Schlcgcl, Finkb eTner fe Davis|199"8| l. We have assumed no 
reddening in the SN host galaxy (see Section 2.3). 



Table 2. PAIRITEL Observations of SN 2008iy. 



'mi 


§ 




J ma 




] 


H ma 






K s 


mag 


h 


(MJ 






(Vega) 






(Vega) 




(V 


ega) 




54934 


42 


16 


59 


± 





05 


16 


00 


± 





05 


14.98 


± 





05 


54937 


43 


16 


57 


± 





06 


15 


99 


± 





06 


15.05 


± 





05 


54941 


43 


16 


57 


± 





05 


16 


04 


± 





08 


15.06 


± 





08 


54944 


40 


16 


70 


± 





05 


15 


98 


± 





07 


14.98 


± 





07 


54948 


10 


16 


59 


± 





05 


15 


87 


± 





06 


14.98 


± 





06 


54951 


39 


16 


55 


± 





04 


15 


96 


± 





06 


15.02 


± 





06 


54954 


35 


16 


70 


± 





04 


15 


88 


± 





05 


15.00 


± 





08 


54957 


37 


16 


63 


± 





05 


15 


91 


± 





05 


14.99 


± 





10 


54959 


37 


16 


68 


± 





05 


15 


89 


± 





06 


14.97 


± 





07 


54961 


11 


16 


69 


± 





06 


15 


99 


± 





08 


14.83 


± 





12 


54970 


30 


16 


77 


± 





08 


15 


98 


± 





12 










54975 


28 


16 


66 


± 





07 


15 


86 


± 





05 


14.87 


± 





06 


54981 


35 


16 


72 


± 





04 


15 


90 


± 





05 


14.93 


± 





08 


54986 


29 


16 


61 


± 





05 


15 


85 


± 





06 


15.06 


± 





12 


54990 


30 


16 


68 


± 





08 


16 


00 


± 





11 


14.59 


± 





07 


54994 


30 


16 


69 


± 





05 


15 


89 


± 





04 


14.84 


± 





06 


54998 


31 


16 


67 


± 





05 


15 


92 


± 





07 


14.90 


± 





07 


55003 


18 


16 


86 


± 





07 


15 


97 


± 





05 


14.96 


± 





11 


55030 


20 


16 


92 


± 





19 


16 


01 


± 





15 










55041 


18 


16 


78 


± 





07 


15 


94 


± 





08 


14.82 


± 





10 


55091 


10 


17 


08 


± 





09 


15 


94 


± 





07 


14.74 


± 





08 



Table 3. Nickel Observations of SN 2008iy. 



date 
(MJD) 


B ma 
(Vcg 


P 


V ma 
(Veg 


a 

g 
a) 


D a 

H mag 
(Vega) 


1 mag 
(Vega) 


54939.52 


19 


06± 





01 


18 


59± 





03 


17 


86± 





01 


17 


87± 


0.03 


54940.46 


19 


09± 





03 


18 


56± 





02 


17 


86± 





02 


17 


89± 


0.03 


54971.37 


19 


35± 





06 


18 


71± 





04 


17 


99± 





01 


17 


98± 


0.11 


54975.36 


19 


25± 





03 


18 


75± 





02 


17 


96± 





02 


17 


92± 


0.07 


54979.46 


19 


23± 





03 


18 


71± 





02 


17 


97± 





02 








54999.35 


19 


36± 





05 


18 


79± 





03 


18 


07± 





03 


18 


12± 


0.08 


55004.37 


19 


40± 





05 


18 


84± 





03 


18 


07± 





02 


18 


13± 


0.11 


55007.43 


19 


56± 





11 


18 


84± 





10 


18 


02± 





01 


18 


37± 


0.13 


55015.36 


19 


43± 





11 


18 


95± 





06 


18 


10± 





04 


18 


26± 


0.08 


55019.30 


19 


58± 





09 


19 


02± 





06 


18 


19± 





03 


18 


21± 


0.07 


55025.36 


19 


66± 





24 


18 


66± 





10 


18 


16± 





08 


18 


18± 


0.12 


55034.26 


19 


52± 





02 


19 


02± 





03 


18 


19± 





02 


18 


23± 


0.11 


55040.30 


19 


56± 





05 


18 


97± 





03 


18 


18± 





03 


18 


30± 


0.09 


55042.24 


19 


60± 





07 


19 


11± 





06 


18 


21± 





03 


18 


31± 


0.11 


55047.27 


19 


43± 





11 


19 


21± 





10 


18 


22± 





01 


18 


24± 


0.20 


55071.22 


19 


66± 





04 


19 


17± 





05 


18 


23± 





03 


18 


34± 


0.09 


55077.23 


20 


07± 





24 


19 


33± 





17 


18 


29± 





03 


18 


58± 


0.17 


55100.15 


19 


75± 





07 


19 


25± 





08 


18 


35± 





03 


18 


49± 


0.06 



a Observed value; not corrected for Galactic extinction. 



a Midpoint between the first and last exposures in a single 
stacked image. 

" Observed value; not corrected for Galactic extinction. 



The reddening of SN 2008iy by the host galaxy is un- 
certain. In our Kast spectra we do not detect any absorption 
from Na I D, and therefore we adopt ^4,\,host = mag. 

We confirm the possible detection of the host galaxy by 
|Catelan et al~| ( [20091 ) in a stack of SDSS g, r, and i-band 
images of the field. When the g + r + i stack is calibrated 
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Figure 2. Filtered photometry including BVRI from the Nickel 
telescope and J, H, and K s from PAIRITEL, showing the decline 
rate of SN 2008iy between days ~560 and 715. To help guide the 
eye, we illustrate the linear fit to each band used to determine the 
photometric decline rates (see text). Note that the single epoch 
of UVOT B and V-band observations we measured, which agree 
to < lcr with Nickel data taken < 2 days later, is not shown in 
this figure. 



Table 4. UVOT Observations of SN 2008iy. 





filter 


ma|] 


^mag 


(MJD) 




(Vega) 




54936.31 


UVW2 


19.14 


0.06 


54936.32 


UVM2 


18.95 


0.08 


54936.32 


UVW1 


18.53 


0.07 


54936.31 


U 


18.60 


0.06 


54937.49 


UVW2 


19.12 


0.04 


54937.85 


B 


19.04 


0.06 


54937.85 


V 


18.65 


0.08 



a Midpoint between the first and last exposures in a single 
stacked image. 

" Observed value: not corrected for Galactic extinction. 



to the SDSS r band, we find that the host has r — 22.75 ± 
0.15 ± 0.08 mag in a 2.2" diameter aperture, where the two 
errors represent the statistical and calibration uncertainties, 
respectively. 



As first noted by Catelan et al. (20091, the field of 



SN 2008iy was observed by the Galaxy Evolution Explorer 



Table 5. Log of spectroscopic observations. 



Epocrj^] 


UT Date 


Telescope / 


Exposure 




Instrument 


« 


560 


2009-04-18.458 


Shane 3-m/Kast 


1500 


567 


2009-04-25.479 


Shane 3-m/Kast 


1500 


652 


2009-07-23.274 


Shane 3-m/Kast 


1800 


711 


2009-09-22.236 


Keck I 10-m/LRIS 


600 



a Defined as rest-frame days relative to day 0, 2007 Sep. 13. 
© 2009 RAS, MNRAS 000,|rjfT5l 



(GALEX; |Martin et al~|p005| . On 2004 May 17, the field 
was imaged as part of the all-sky imaging survey (AIS); 
the host was not detected to a 5a limiting magnitude of 
FUV > 20.2 and NUV > 20.5 mag. The field was reimaged 
as part of the medium imaging survey on 2008 Jun. 05/06, 
and SN 2008iy was detected at FUV = 21.175 ± 0.094 mag 
and NUV = 19.929 ± 0.036 mag. 



3 RESULTS 

3.1 Photometric Analysis 

At z = 0.0411, the absolute peak magnitude of SN 2008iy 
was a relatively modest M r « —19.1 mag. This places 
SN 2008iy well below the most luminous SNe Iln, such as 
SN 2006gy (Mr = -21.7 mag; |Ofek et al."|[2007 



et al. 



et al. 



20071 and SN 2008fz (M v 



-22.3 mag; 



Smith 



Drake 



2009b I, and more in the range of the well-studied 



SN Hn 1988Z (Mr < -18.9 mag; Turatto et al. 1993 



Stathakis & Sadler||1991|). Assuming no bolometric correc- 



tion, the total integrated optical output from SN 2008iy dur- 
ing the first ~700 days after discovery is ~2 x 10 50 erg. Be- 
tween day 550 and 720 the SN declines in the optical and 
the J band, while it actually gets brighter in the H and 
K s bands. During this time, the linear decay rates are as 
follows: /3b = 0.51 mag (100 day)" 1 , j3 V = 0.50 mag (100 
day)" 1 , Pr = 0.34 mag (100 day)" 1 , /3 7 = 0.43 mag (100 
day)" 1 , /3j = 0.30 mag (100 day)" 1 , j3 H = -0.03 mag (100 
day) -1 , and /3k s = —0.19 mag (100 day) -1 . These decline 
rates, all slower than the expected rate of decline for ra- 
dioactive 56 Co, 0.98 mag (100 day) -1 , strongly suggest that 
the SN is still being powered by CSM interaction ~700 days 
after explosion. 

For comparison purposes, in addition to the light curve 
of SN 2008iy, Fig. [I] s hows the light curve of the standard 
Type II-P SN 1999em ( |Leonard et al. |2002[ > as it would have 
appeared at the redshift of SN 2008iy. As well as illustrating 
the very long rise time of SN 2008iy, this comparison shows 
that SN 2008iy was fairly luminous at the time of discov- 
ery. Thus, despite the sparse sampling over the first ~200 
days, this large luminosity suggests that the early detections 
are not related to a pre-SN outburst, as was observed for 
SN 2006jc ( |Foley et al. [2007| |Pastorello et al. |200T| ). CRTS 
images indicate that the SN was rising over the first ~400 
days ( Mahabal et al. |2009 1 , which provides further evidence 
against a pre-SN outburst. Fig. [j] also shows the light curve 
of the long l ived, interacting, Type Iln SN 1988Z ( |Turatto| 



et al. 



1993[ ). The explosion date of SN 1988Z is not well 
constrained (see Section 4.3), so we shift the first epoch of 
detection to day 0. Notice the similarity of the decline rate 
of SN 2008iy and SN 1988Z around day ~600; these SNe 
also have very similar late-time spectra (see Section 4.2). 

We were unable to fit the UV through NIR spec- 
tral energy distribution (SED) of SN 2008iy to a single- 
temperature blackbody model. This is not surprising, as 
SNe II are dominated by emission lines at late times, and 
single-temperature blackbody models typically apply only to 
young SNe II (see, e.g., Filippenko 1997). Direct integration 
of the UV— NIR SED shows that the bolometric correction, 
relative to the I band, is a factor of ~3 in luminosity, cor- 
responding to ~1.2 mag, on day ~560. Some features stand 
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out from the SED: there is a strong i?-band excess rela- 
tive to the other optical bands, with V — R = 0.8 mag and 
R — I = —0.2 mag on day 715. The red V — R colour relative 
to a blue R — I colour can be attributed to the Ha emission 
with large equivalent width. This emission also accounts for 
the relatively slow decay rate in the R band, as compared to 
the other optical bands. There is also a NIR excess relative 
to the J-band flux. In fact, this excess increases with time, 
I — K m 2.8 mag on day ~560 and / — K « 3.7 mag on day 
~710, which is indicative of the growing importance of dust 
in the emission from SN 2008iy. 

A NIR excess at late times has been observed in many 
SNe Iln (e.g., Gerardy et al. |[2002 |, and it is the result of 
either new dust formation (e.g., SN 2005ip; Smith et al 
2009| |Fox et al. ||2009[ ) or the pres ence of a NIR light echo 
from preexisting dust (Dwck 1983), or both. To distinguish 
between these two possibilities, which are virtually identi- 
cal from photometry alone, requires well-sampled optical 
spectra (see Section 3.2.2) since line profiles are expected 
to change with time if new dust is being formed. On day 
706, corresponding to our last PAIRITEL observation, the 
H — K s colour of SN 2008iy was ~1.2 mag. Assuming that 
the dust radiates as a perfect blackbody, this colour corre- 
sponds to a dust temperature of Td ust ~ 1320 K, while 
the i^-band measurement corresponds to a dust luminosity 
Tvnir ~ 4.8 x 10 42 erg s _1 , assuming no bolometric cor- 
rection as we cannot constrain the emission in the mid-IR. 
This luminosity is very large, though upon making similar 



assumptions Gerardy et al. (2002) found NIR luminosities 
> 10 42 erg s 1 for SNe 1995N and 1997ab at late times. Ger- 
ardy et al. also found that the NIR luminosity was roughly 
an order of magnitude greater than the X-ray emission from 
SN 1995N, which is also the case for SN 2008iy (see Sec- 
tion 3.3). The H — K s colour of SN 2008iy is bluer than 
those of the Gerardy et al. sample at a similar epoch, which 
may be a result of the long rise of SN 2008iy or a contri- 
bution to the NIR emission from the SN in addition to the 
dust. Future NIR observations, as the underlying SN light 
continues to fade, will place more stringent constraints on 
the dust near SN 2008iy. 

3.2 Spectroscopic Analysis 

Our spectra of SN 2008iy at ages >560 days are very simi- 
lar to those of an unpublished spectrum, from our spectral 
database, of the Type Iln SN 1988Z taken at a compara- 
ble epoch, as shown in Fig. [3] The general features resem- 
ble those of several other SNe Iln ( |Filippenko|| 1997[ ) : there 
are prominent Balmer and He I emission lines, most no- 
tably the dominant Ha emission, which primarily feature 
intermediate-width emission components with full width at 
half-maximum intensity (FWHM) ~1650 km s" 1 in the case 
of SN 2008iy. The higher-resolution Keck spectrum reveals 
a number of narrow, marginally resolved (FWHM < 170 km 
s" 1 ) emission lines, including Ha, H/3, [O III] AA4363, 4959, 
5007, and He I A7065. The only high-ionisation line we detect 
is [Fe X] A6375. The relative lack of narrow forbidden lines 
and the low intensity ratio of [O III] AA4959, 5007 to [O HI] 
A4363 suggest a large electron density f or the ejecta (|Filip- 



We can estimate the density of the unshocked emit- 
ting material based on the relative intensities of the three 
[O III] lines mentioned above. Note that [O III] A4363 is 
only resolved in our day 711 Keck spectrum, so the follow- 
ing estimate of the density is at that epoch only. Following a 
removal of the underlying continuum, we fit three Gaussian 
profiles to [O III] A4363 and [O III] AA4959, 5007, and find 
that TZ = /[A4959 + A5007]/J[A4363] « 1.7. As mentioned 
above, the strong presence of [O III] A4363 indicates a high 
density for the emitting material. In fact, with TZ ~ 1.7, the 
electron density, n c , must be > 10 6 cm -3 regardless of the 
temperature of the emitting gas. Typically, [O III] emission 
comes from photoionised regions with T = 1 6, 000-20,000, 
in which case n e « 10 7 cm -3 (see Fig. 11 in Filippenko & 



Halpern|1984 1. Note that this high-density material is likely 



only present within clumps in the CSM (see Section 4.2). 
and that interclump portions of the CSM have a lower den- 
sity. 

3.2.1 The Ha Profile 

As seen in Fig. [3j the Ha emission feature dominates over all 
the other lines. In the high-resolution Keck spectrum, we see 
evidence for three distinct emission features: a broad compo- 
nent (FWHM « 4500 km s" 1 ), an intermediate component 
(FWHM « 1650 km s _1 ), and a marginally resolved, narrow 
P Cygni feature (FWHM » 75 km s" 1 ), as shown in Fig. [Z] 
The main panel of Fig. [4] shows a fit to the Ha profile (red 
line) which includes a broad (FWHM w 4500 km s _1 ) and an 
intermediate (FWHM « 1650 km s -1 ) component. Panel (a) 
in Fig. [4] shows a close-up view of the narrow emission after 
the broad and intermediate components have been removed 
using a spline fit. We fit the spline to the observed Ha profile 
from —1000 km s -1 to 1000 km s -1 after masking the region 
between —200 km s" 1 and 200 km s _1 . A narrow absorp- 
tion minimum such as this, located at ~ —100 km s~ , has 
been seen in a number of typical SNe Iln (e.g., SNe 1997ab, 
1997eg, and 1998S, |Salamanca, Terlevich fe Tenorio-Taglel 



penko fe Halpern| 1984| [Filippenko| 1989) |Stathakis & Sadler 



1991 ). The relative spectral features do not show strong evo- 



lution between days 560 and 711. 



2002||Salamanca|2000||Fassia et al. |2001[ ), and it traces the 
outflow velocity of the progenitor wind. We note that while 
the narrow absorption is only marginally resolved, a P Cygni 
profile with characteristic speeds of ~10 km s" 1 , typical of 
red supergiants (RSGs), would go completely unresolved in 
the Keck spectrum. Therefore, this P Cygni profile must 
be associated with an outflow velocity that is >10 km s~\ 
which will have important consequences for the progenitor 
(see Section 4.3). 

It is interesting that the broad plus intermediate fit 
overestimates the Ha flux out to about —450 km s~ . In 
panel (b) of Fig. [4] we show the residual flux, centred on 
Ha, following the subtraction of our two-Gaussian fit to the 
broad plus intermediate emission. The blue velocity at zero 
intensity (BVZI) extends to roughly —450 km s _1 , which 
means that the true wind speed from the progenitor, as 
traced by the absorbing gas moving directly along the line 
of sight toward the SN, may be >100 km s _1 and possibly 
as high as ~450 km s _1 . If this feature represents a true 
lack of emission, it would be evidence for a two-component 
CSM: one with velocity u w ,i ~ 100 km s _1 and the other 
with BVZI w n w ,2 450 km s~ . Similar features were seen 
in SN 1998S (|Fassia et al. ||200l| , which were modelled by 
Chugai et al. ( 2002 1 to be a slow wind accelerated to higher 
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Figure 4. Detailed view of the Ha profile from the Keck spectrum 
taken on day 711. Main panel: The black line shows the observed 
Ha profile, while the smooth red line illustrates a two-Gaussian 
fit (intermediate, FWHM « 1650 km s _1 ; broad, FWHM « 4500 
km s _1 ) to the profile. Residual emission and absorption from a 
narrow component to the profile can easily be seen. The location 
of He f A6678 has been marked, and the residuals from an im- 
perfect telluric absorption correction have been marked by the © 
symbol. Inset (a): Close-up view of the residual, narrow P Cygni 
component. The smooth red line shows a fit to the feature that in- 
cludes a Gaussian in emission and a Gaussian in absorption (both 
Gaussians have FWHM PS 75 km s — x ). The absorption minimum 
of the narrow absorption is located at —100 km s _1 . Note that 
for this fit we removed the broad plus intermediate components 
using a spline fit over the region from —1000 km s~ 1 to 1000 
km s —1 (see text). Inset (b): Residuals following the subtraction 
of the two-Gaussian broad plus intermediate fit from the observed 
Ha profile. The absorption minimum still occurs at —100 km s _1 , 
but the BVZI extends to ~ —450 km s _1 , meaning the progenitor 
wind speed may be > 100 km s" 1 (see text). 



velocities by radiation from the SN photosphere, and in SN 
2006gy, which Smith et al. ( 2010 1 argued was the result 



of a CSM shell that had been ejected from the progenitor 
in a pre-SN eruption. These two scenarios predict different 
behaviour as the SN evolves. 

In the radiatively accelerated scenario, as shown by 
|Chugaiet al~| ( |2002[ ) for SN 1998S, the second, faster compo- 
nent to the wind has a negative velocity gradient. Observa- 
tionally, this effect manifests itself as a BVZI for the second, 
faster CSM component that decreases with time, as observed 
in SN 1998S ( |Fassia et al~||2001| . In the shell-ejection sce- 
nario, the wind exhibits a positive velocity gradient as the 
ejected shell is freely expanding, which is thought to be the 
result of an explosive (~ 10 49 — 10 50 erg) mass-loss event (see 
e.g., |Chugai et al. |2004| ). This scenario leads to a BVZI that 
increases with time, as was observed for SN 2006gy (Smith 



et al. 20101. With only a single high-resolution spectrum, 



we are unable to probe the evolution of this feature, but we 
note that the veloc ity of a radiatively accelerated CSM i s 



proportional to t (Fransson, Lundqvist & Chevalier 1996), 



so this mechanism is unlikely to be significant at late times 
(> a few hundred days). 

As does the overall continuum, the Ha luminosity de- 
clines over the course of our observations. In Table [6] we 
summarise the observed properties of the Ha profile from 
our spectra on days 560, 652, and 711. We do not include 



Figure 5. Evolution of the Ha profile of SN 2008iy between day 
560 and day 711. The profile remains largely unchanged, with the 
exception of a variable spectral feature between ~200 km s — 1 and 
20,000 km s — 1 , from day 652 to day 711. Similar spectral features 
have been observed in other interacting SNe (see text). Emission 
lines of He I A6678 and He I A7065 are labelled. 



the spectrum from day 567, as it was taken under cloudy 
conditions and the overall flux scaling is uncertain, but we 
do not expect significant spectral evolution between days 
560 and 567. The decline in the Ha luminosity is accompa- 
nied by a rise in the equivalent width over this same time 
period. For the narrow emission on day 711, we measure 
L(Ha) = (2.8 ± 0.6) x 10 39 erg s~\ though we note that 
this value is likely underestimated because the blue emission 
wing is probably being partially absorbed. The FWHM of 
the intermediate component is ~1650 km s" 1 at all epochs, 
while the FWHM of the broad component drops from ~5800 
km s _1 on day 560, to ~5200 km s" 1 on day 650, to ~4400 
km s" 1 on day 711. We attribute this change in the broad 
component to rapid evolution of the spectrum redward of 
Ha between days 560 and 711, as shown in Fig. [5] The na- 
ture of this feature is not currently understood: it could 
potentially be an absorption feature, or it could be related 
to a lack of emission from some unidentified species. Similar 
spectral evolution redward of Ha has been seen in the inter- 
acting SN e 2005gj and 2006gy ( |Prieto et al. |2007[ [Smlthetl 
al. |2010 l. We do not associate this feature with Ha, as the 
red velocity at zero intensity corresponds to velocities of > 



22,000 km s , and there are no other emission lines that 
show evidence for velocities this large. 



3.2.2 He I Emission 

In addition to the prominent Balmer emission lines, 
SN 2008iy exhibits several He I lines in emission, includ- 
ing AA3820, 4026, 4471, 5016, 5875, 6678, and 7065. None 
of these lines show evidence for a broad component, and the 
intermediate component, FWHM « 1100-1400 km s -1 , is 
slightly narrower than that observed for Ha. As previously 
mentioned, He I A7065 has a narrow, marginally resolved 
emission component in addition to the intermediate compo- 
nent. 

A systematic, increasing blueshift in the He I profiles 
has been argued as evidence for dust formation in the cool, 
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Table 6. Ho line widths, luminosity, equivalent width, and ratio to H/3. 



EpoclQ 


Br. FWHM 


Int. FWHM 


Nar. FWHM 


L(Raf\ 




Ha/H/3 


(km s —1 ) 


(km s _1 ) 


(km s -1 ) 


(erg s ^) 






560 


5774 


1641 




7.86 x 10 41 


2122 


9.49 


652 


5162 


1685 




6.10 x 10 41 


2386 


12.22 


711 


4405 


1639 


78 


6.50 x 10 41 


2786 


12.54 



a Defined as rest-frame days relative to day 0, 2007 Sep. 13. 
b Uncertainty ~ ±10%. 
c Uncertainty ~ ±10%. 



dense shell formed in the post-shock gas of the Type Iln SNe 
1998S and 2005ip ( [Pozzo et al~|[2004l [Smith et al. p009 |. 
Newly formed dust in the post-shock region absorbs light 
from the receding ejecta, resulting in a suppression of flux 
from the red side of the emission line. Using He I A7065, the 
strongest He line in our spectra, we searched for changes 
in the red wing of the profile. The Kast spectra from days 
560 and 567 are very noisy, but in the higher signal-to-noise 
ratio spectra from days 652 and 711 we see no evidence for 
a change in the red wing of the emission profile. This hints 
that the observed NIR excess (see Section 3.1) is the result 
of a NIR echo, rather than newly formed dust, though we 
caution that with only two spectra separated by ~50 days 
the lack of change in the He I profile does not constitute 
definitive evidence for this. 



3.3 X-ray Analysis 

To determine the X-ray luminosity of SN 2008iy we assume 
a thermal plasma spectrum with temperature 0.6 keV (see, 

which has been abso rbed by a 
2 dKalberla 



e.g 



Immlcr & Kuntz 2005) 



Galactic column density of n(H) = 4.8 x 10 20 cm 
|et al. |[20fJ5) . Using the 0.3-10 keV XRT count rate (see 
Section 2.1) and PIMM^J we measure the unabsorbed flux, 
assuming no absorption within the host, to be (6.1 ± 1.9) 
x 10~ 14 erg cm~ 2 s" 1 . At the adopted distance of SN 2008iy 
this corresponds to a luminosi t y of L x = (2.4 ± 0.8) x 10 41 
erg s" 1 . Note that |Fox et al. | ( 2000| fit the late-time (> 3 



yr) X-ray spectra of SN 1995N with a larger temperature of 
9.1 keV, which would correspond to a flux and luminosity of 
(9.1 ± 2.9) xlO" 14 erg cm -2 s _1 and (3.5 ± 1.1) x 10 41 erg 
s _1 , respectively, were we to adopt this value for SN 2008iy. 
Similar luminosities (Lx > 10 41 erg s _1 ) are observed for 
other interacting SNe more than 1 year post explosion (e.g., 
SNe 1988Z and 1995N; |Fox~et al. |2000[ ). 



4 DISCUSSION 

4.1 Mass Loss in the SN 2008iy Progenitor 

We assume that the continuum emission traced by the light 
curve of SN 2008iy is powered solely by CSM interaction. 
The half life of 56 Ni is r ~ 6.0 days, much too short to 
power the > 400 day rise of SN 2008iy, while the photometric 
decay rates are too slow to be powered by 56 Co (see Sec. 

5 http:/ /heasarc. nasa.gov/Tools/w3pimms. html. 



3.1), implying that radioactivity contributes minimally to 
the light curve of SN 2008iy. Furthermore, the SN ejecta 
will suffer considerable adiabatic losses during the 400 day 
rise of the light curve, suggesting that interaction is the only 
viable mechanism for powering SN 2008iy over this entire 
time. 

Following our assumption of interaction-driven lumi- 
nosity, the expelled gas from the progenitor overtaken by 



the shock at each epoch can be written as (see Chugai & 
|Danziger||l994| ) 



(1) 



SN 



where M is the mass-loss rate of the progenitor, ij) is an 
efficiency factor describing the conversion of kinetic energy 
into radiation, L is the luminosity of the SN, v v is the wind 
speed of the progenitor, and Vsn is the velocity of the blast 
wave overrunning the CSM. We adopt ij) = 0.5, but this 
number is likely an overestimate given the optically thin na- 
ture of the emission (see Section 4.2). The CSM wind is only 
probed by our higher-resolution spectrum from day 711, so 
we adopt « w = 100 km s" 1 based on the narrow P Cygni 
absorption seen in Ha, though this may not be valid at all 
times between days and 700 (see below). We adopt Vsn = 
5000 km s _1 based on the typical widths of the broad Ha 
emission, which should trace the speed of the blast wave. 
Unlike our assumption of a constant wind speed, the adop- 
tion of a constant-velocity blast wave at all times should at 
the very least be valid after the peak, at which point the 
CSM density is decreasing with increasing radius. Hydro- 
dynamic modelling of SNe Iln demonstrates that the blast 
wave is quickly decelerated to a roughly constant expansion 
speed as the ejecta sweep up successively more CSM (see, 
e.g., Chugai et al. 20041. This scenario may not apply to 
SN 2008iy, however, because the ~400 day rise time implies 
that the CSM density may have been increasing with radius, 
in which case the blast wave may have been continually de- 
celerated during the rise. 

Following the above assumptions we can use Equation[T] 
to determine the mass loss traced by the optical luminosity 
at a number of different epochs of interest. At each epoch 
we determine the radius, R = Vsn£sn ~ 5000 km s" 1 tsn, 
and the luminosity based on the light curve shown in Fig. [I] 
We do not adopt the bolometric correction determined in 
Section 3.1 for two reasons: (i) it is unclear whether this 
same correction is valid at early times, and (ii) emission from 
a NIR echo should not be incorporated into Equation [T] as 
that contribution to the total luminosity is not directly the 
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result of CSM interaction. Thus, the estimates below of the 
luminosity do not reflect any emission from dust. At early 
times, day ~153, L « 6.8 x 10 8 Lq, which corresponds to 
M « 1.3 X 10~ 2 M yr" 1 at R « 6.6 x 10 15 cm, while at 
the time of peak L w 1.2 x 10 9 Lq, R « 1.7 X 10 16 cm, and 
M w 2.3 x 10~ 2 M yr -1 . On day 563, coincident with our 
X-ray observations of SN 2008iy the optical luminosity was 
isi7x 10 8 L Q , which corresponds to R « 2.4 x 10 16 cm and 
M « 1.3 X 10" 2 M Q yr -1 . On day 714, roughly coincident 
with our high-resolution Keck spectrum, L w 4 x 10 s L©, 
7? « 3.0 x 10 16 cm, and M « 0.8 x 10~ 2 M Q yr" 1 . 

These results are quite remarkable; they suggest that 
at a time t « Vsn^sn/uw « 5000 km s" 1 x 400 day/100 
km s _1 w 55 years prior to the explosion of SN 2008iy, the 
progenitor underwent a period of heightened mass loss. This 
enhanced mass-loss period was then followed by a period of 
decreasing mass loss leading up to the time of the SN ex- 
plosion. While the quantitative results for the mass loss dis- 
cussed above are sensitive to our adopted quantities, most 
specifically the progenitor wind speed and the speed of the 
SN blast wave, the fact remains that the continuum lumi- 
nosity did increase over a period of ~400 days after the 
SN explosion. Rearrangement of Equation [T] shows that this 
luminosity increase is proportional to the wind-density pa- 
rameter, id = M/w w . Thus, w must have increased over a 
distance R ~ Vsn^sn ~ 1.7 x 10 16 cm from the progenitor, 
regardless of the true wind speed over that distance. An in- 
creasing value of w means that during roughly the century 
prior to the SN, there was a significant change in the wind 
properties of the progenitor. 

Using our X-ray detection of SN 2008iy we have an al- 
ternative method to probe the mass-loss history of the pro- 
genitor on day ~563. Following Immler & Kuntz ( 2005 1, the 



X-ray luminosity may be written as 



M 



{■Km)' 2 



A(T) - (Vsnty 1 , 



(2) 



where Lx is the X-ray luminosity, m is the mean mass per 
particle (m = 2.1 x 10~ 24 g for a H+He plasma), A(T) is the 
cooling function of a plasma heated to temperature T, M is 
the mass-loss rate of the progenitor, v v is the wind speed of 
the progenitor (~ 100 km s _1 , see above), Vsn is the speed 
of the SN blast wave (~ 5000 km s _1 , see above), and t is 
the time since explosion, ~563 days. Assuming a tempera- 
ture T = 10 7 K, appropriate for an optically thin thermal 



plasma (Immler & Kuntz 20051, the effective cooling func- 



tion is A(T) = 3 x 10 _: " erg cm 3 s . The X-ray luminosity 
on day 563 was Lx ~ 2.4 x 10 41 erg s _1 . Substituting these 
values into Equation^ we find M » 0.7 x 10" 2 M yr" 1 , 
which shows agreement with the value derived from the opti- 
cal continuum, ~ 1.3 x 10 -2 M© yr -1 , given the uncertainty 
in a number of the assumptions we have adopted. 

An additional estimate of the mass-loss rate of the pro- 
genitor comes from the narrow Ha? emission seen in our Keck 



spectrum on day 711. Chugai & Danziger (2003) show that 



the Ha emission from the unshocked CSM is related to w. 

L(Ha) = — — as2 hu 23 {xXwN A ) 2 f 1 - — ) , (3) 

where ri is the inner radius corresponding to the position of 
the SN blast wave, r 2 is the outer radius related to the fast- 
moving forward shock, 0132 is the effective recombination 



coefficient for Ha, hv 23 is the energy of an Ha photon, x is 
the degree of H ionisation, X is the hydrogen mass fraction, 
and Na is Avogadro's number. We cannot constrain x, X, or 
r 2 , so we conservatively adopt x = 1, X = 1, and r 2 3> n, 
all of which have the effect of minimising w. Assuming Case 
B Ha recombination, 032 = 8.64 x 10~ 14 cm 3 s~\ which is 
appropriate assuming that the narrow emission comes from 
photoionised gas. Substituting the narrow Ha luminosity 
(L(Ha) = 2.8 x 10 39 erg s" 1 ) and n = Vsntsn = 3.0 x 10 16 



cm into Equation I3I we find that w w 10 1 ' g cm -1 . With 
v w — 100 km s~ on day 711, this corresponds to M ~ 
1.6 x 10~ 2 M© yr -1 . At a similar epoch (see above), we find 
M « 0.8 x 10 -2 M© yr -1 based on the optical luminosity. 
Again, the agreement between these two methods to within 
a factor of ~2 is reasonable given the uncertainties in our 
adopted parameters. 



4.2 Late-Time Emission and the Similarity to SN 
1988Z 

To explain the late-time (i.e., post peak) emission from 
SN 2008iy we adopt a model that is virtually identical to 



that developed by Chugai & Danziger ( 1994 1 for SN 1988Z 



the CSM contains optically thick clumps in addition to a rar- 
efied wind between these clumps (see, e.g., Fig. 1 of Chugai 
fc Danziger|1994| >. As the SN ejecta sweep through the CSM, 



they drive a fast-moving forward shock into the rarefied 
wind. This shocked material forms a cool dense shell between 
the forward and reverse shocks, and gives rise to the broad 
emission component seen in Ha (FWHM w 5000 km s _1 for 
SN 2008iy) . At the same time a slower shock is being driven 
into the dense clumps, which leads to the intermediate- width 
emission (FWHM w 1650 km s" 1 in the case of SN 2008iy) 
seen in the spectra. The narrow P Cygni profile results from 
the pre-shock photoionised wind. 

The observed similarities between SN 2008iy and 
SN 1988Z provide further evidence that comparable mod- 
els are appropriate for the two SNe. In addition to hav- 
ing similar spectra at late times (see Fig. |3|, SNe 1988Z 
and 2008iy have large X-ray luminosities (both SNe have 
Fab>ian fc Terlevich|l996[ this work), and 



10 41 ergs- 



similar late-time decline rates (~0.3— 0.5 mag (100 day) -1 
in the optical; [Turatto et al. |1993| this work), which in both 
cases are slower than the expected bolometric decline of ra- 
dioactive 56 Co. Estimates for the mass-loss rate from the 
progenitor of SN 1988Z vary by roughly an order of mag- 
nitude. From X-ray observations Schlegcl & Pctre 
find A/1988Z w 10~ 3 M© yr" 1 , while |Williams et al. 
estimate Mi9gsz ~ 10 -4 M© yr -1 based on radio observa- 
tions, and the m odels of the optical emission by Chugai & 



2006 



2002 



Danziger ( 1994 1 yield Mi 988 z « 7 x 10" 4 M© yr -1 . While 



these estimates are 1-2 orders of magnitude less than those 
for SN 2008iy, we note that in each of the above estimates 
for Migssz a wind speed of 10 km s _1 was adopted for the 
progenitor of SN 1988Z. If SN 1988Z had a progenitor wind 
speed closer to ~100 km s _1 , a very reasonable possibility 
given the wind speed of SN 2008iy and other similarities be- 
tween the two SNe, this would result in an increase in the 
estimated mass-loss for the progenitor by a factor of ~10, 
thereby bringing the estimates for the two SNe into accord. 

Based on the late-time emission, SN 2008iy seems to 
belong to the group of SNe Iln that exhibit a slow evolution 
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sustained by long-lived (> 1 decade) CSM interaction. In 
addition to SN 1988Z, other members of this group include 
SN 1986 J (jRupen et al. (1987] ), SN 1995N ( |Fox et al. |2000 



Fransson et al. |2002| >, and possibly also the VLSN 2003ma 



(Rest et al. 20091. Both SN 1986 J and SN 1995N exhibit 



evidence for a clumpy progenitor wind, like SN 2008iy and 



SN 1988Z. Chugai (19931 modelled the X-ray emission from 



century prior to explosion, this could result in a density pro- 
file that peaks ~1.7 xlO 16 cm from the progenitor. In the 
clumpy-wind scenario described in Section 4.2, this would 
mean that the ejecta are expanding into a wind where the 
number density of clumps is increasing with radius. Thus, as 
more and more clumps are overtaken by the ejecta, the con- 
tinuum luminosity continues to rise, until after ~400 days 
when the ejecta have reached ~1.7 xlO 16 cm, the number 
density of clumps begins to decline and so does the optical 
luminosity. Alternatively, the long rise could result from the 
ejecta expanding into a nonspherical wind, such as a bipolar 
outflow, though this hypothesis would need to be examined 
with detailed hydrodynamical models. 

The inferred mass-loss rate for SN 2008iy, M « 10 -2 
Mq yr -1 , is similar to that for the first LBV, P Cygni, dur- 
1 mag in the optical between 1994 and 2003 ( Milisavljevic et ing its great eruption (Mpc yg ~ 10 -2 Mq yr -1 during the 



SN 1986J as the interaction between the SN ejecta and a 
clumpy wind (though note that Houck et al. |[l998| prefer a 
model with a smooth CSM, but they conclude that they can- 
not rule out the clumpy model). The evidence for a clumpy 
wind from the progenitor of SN 1995N comes from both the 
optical spectra (Fransson et al. | 2002[ ) and the X-ray emis- 
sion ( Zampieri et al. |2005 1. Also, the optical decline of both 
SNe 1986J and 1995N is very slow: SN 1986 J declined by < 



|al. ||2008| >, while the V-band dec line of SN 1995N was only 
2.2 mag between 1998 and 2003 ( |Zampieri et~al~|[2005 1 . 

The duration of the interaction means that the large 
mass-loss rates from the respective progenitors must have 
been sustained for at least ~ 100 yr prior to core collapse 
after the conservative assumption that Vsn = 1000 km s~\ 
and v w = 100 km s _1 . Typically the blast wave continues to 
expand at > 1000 km s _1 (SN 1988Z had a broad component 
that remained nearly constant at ~2000 km s _1 from day 

19991), which would make 



1500 to 3000; Aretxaga et al. 



this time period > 100 yr. These very long-lived SNe Iln are 
therefore connected in that their progenitors experienced 
lengthy periods with high mass-loss rates. This stands in 
stark contrast to the Type Iln SN 1994W, which had a light 
curve with an abrupt drop ~100 days post explosion ( Soller- 
man, Cumming fc Lundqvist[T9 98). The mass loss from the 
progenitor of SN 1994W has been modelled to occur in a 
short (<1 yr), violent episode (Chugai et al. 2004), and 



the sudden drop in the light curve occurs because the SN 
ejecta have overtaken the dense CSM, thereby halting any 
interaction luminosity. 



4.3 Origin of the 400 Day Rise Time and 
Implications for the Progenitor 

With a model to account for the late-time emission from 
SN 2008iy, we are still left with the puzzle of explaining the 
400 day rise time|^] This rise in the optical is significantly 
longer than that seen in any other known SN. As previously 
discussed, this scenario is possible if the progenitor under- 
went a phase of enhanced mass loss ~55 years prior to the 
SN explosion. The optical luminosity traces the wind-density 
parameter, w, meaning that during the decades prior to ex- 
plosion either (i) the mass-loss rate declined, or (ii) the wind 
speed increased, or (iii) both. 

Episodic periods of enhanced mass loss, sometimes in 
the form of a shell ejection, have been observed for numer- 



ous LBVs (see Humphreys & Davidson 19941, though the 
underlying physics of these eruptions is not currently well 
understood (see |Smith fc Owocki|2006[ ) . If the progenitor of 
SN 2008iy had a giant eruption (similar to that of LBVs) < 1 



1600 AD eruption; Smith & Hartigan 2006). Furthermore, 
the terminal wind speed in the P Cygni nebula is 185 km s _1 



( Lamers et al. |1996|[Najarro, Hillier fc Stahl|1997[ ), which is 
similar to the observed BVZI for the narrow blue absorption 
seen in SN 2008iy, corresponding to n w ~ 160-450 km s . 
We illustrate these comparisons to show that derived prop- 
erties of the progenitor wind of SN 2008iy are similar to the 
giant eruption of a Galactic LBV; we are not insinuating 
that the great outburst from P Cygni is a direct analogue 
to the proposed LBV-like eruption from SN 2008iy. Such 
eruptions are not expected from RSGs. 

Further evidence that the progenitor of SN 2008iy could 
not have had a RSG progenitor comes from the observed nar- 
row P Cygni Ha profile. RSGs have typical wind speeds of 
~10 km s -1 , with extrem e RSG winds reaching 40 km s" 1 
(see |Smith et al. |[2007| . The observed 100 km s" 1 wind 



speed from SN 2008iy is more characteristic of LBVs (e.g., 
P Cygni, see above) or the escape velocity from a blue su- 
pergiant. Similar ~100 km s" 1 P Cygni profiles have been 
seen in a numb er of SNe Iln, such as SNe 1997ab ([Salamanca| 
20001), 1997eg (ISalamanca, Terlevich & Tenorio-Tagle|2002l), 



1998S ( |Fassia et al. |2001| ), and 2007rt ( jTrundle et al. |2009| , 
implying that the progenitors of each of these SNe may have 
been in a similar state shortly before core collapse. 

The dense CSM and large luminosities associated with 
many SNe Iln have led a number of authors to suggest 
that at least some SNe Iln are associated with progenitors 
that experienced LBV-like mass loss shortly before core col- 
lapse (see e .g.,|Chu et al. ||1999| |SaIamanca||2000| |Chugai 
fc Danziger]|2003| |Chugai et al. ||2004[). This possible con- 



nection was considerably strengthened following the direct 
identification of the progenitor of the normal Type Iln SN 



2005gl on archival Hubble Space Telescope images ( Gal- Yam 
et aL |2007| |Gal-Yam fc Leonard|2009[ ). While the diagnos- 



tics above are not ubiquitous for all SNe Hij^] many of the 
signatures (high-density CSM, episodic mass loss with M 
~ 10~ 2 Mq yr -1 , and ~100 km s _1 progenitor wind speed) 
are shared with SN 2008iy, suggesting that it too had an 
LBV-like progenitor. 

As mentioned above, an alternative way to generate 
a wind-density parameter that is increasing with distance 



6 We have no spectra of the SN while it was still on the rise, 
thus we have no way of knowing if the SN underwent significant 
spectral evolution prior to day 560. 



7 Note that the low-resolution spectrographs typically employed 
to observe SNe lack sufficient resolution to resolve narrow 
(FWHM < 150 km s _1 ) absorption lines. 
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from the progenitor is to increase the progenitor wind speed 
in the years leading up to the SN event. A relatively mild 
change in the wind speed by a factor of ~3 could, in fact, 
explain the observed increase in luminosity. During the post- 
main-sequence evolution of very massive stars a slightly 
more extreme transition is expected as stars exit the RSG 
or LBV phase, with wind speeds of ~10— 100 km s _1 , to be- 
come Wolf-Rayet stars, with typical wind speeds of > 1000 
km s _1 . As this fast wind expands into the slower wind it 
will sweep up a thin, dense shell at the interface of the two 



winds and thus create a wind-blown bubble (e.g., Dwarkadas 



2005). A SN exploding in one of these bubbles will appear 



as a normal SN until the ejecta reach the edge of the bubble. 
As the ejecta overtake the thin shell the appearance of the 
SN will dramatically change as CSM interaction begins to 
provide a significant contribution to the broad-band lumi- 
nosity. This is the precise scenario proposed for the Type Iln 
SN 1996cr, for which Bauer et al. | ( |2008[ ) found an optical 
luminosity that is clearly decreasing between 1996 and 1999, 
but the radio emission sharply increases in late 1997 while 
X-ray emission abruptly turns on some time between 1998 
Mar. and 2000 Jan. The X-ray emission is especially inter- 
esting because after the abrupt initial rise, it continues to 
rise for at least the next seven years, a feature which had not 
previously been observed for any X-ray SN. A similar sce- 
nario is proposed for SN 2001em, classified as a Type Ic, but 
which showed unexpectedly strong radio and X-ray emission 
> 2 yr after explosion, as well as a FWHM w 1800 km s _1 



Ha emission line (Soderberg, Gal- Yam & Kulkarni 2004) 



Each of these latc-timc peculiarities can be understood if 
the progenitor of SN 2001em lost the remainder of its hy- 
drogen envelope during a phase of large mass-loss ~(2-10) 
xlO" 3 M Q yr" 1 , about (1 -2) xlO 3 yr prior to core collapse 
(Chugai & Chevalier 2006). After the hydrogen envelope has 
been lost, a fast wind from the progenitor sweeps the hydro- 
gen into a shell. The interaction of the Type Ic ejecta with 
this swept-up shell give rise to the unusual late-time emis- 
sion. 

Can the long optical rise be the result of SN 2008iy ex- 
ploding inside a wind-blown bubble? There is a somewhat 
sharp rise in the optical (> 2 mag brightening during the 68 
days between the last DS non-detection and the first CRTS 
detection of the SN) followed by a slow and steady increase 
over the next ~400 days, which may hint that the SN ex- 
ploded in such a wind-driven cavity. However, we do not 
favour this model. In the wind-blown bubble scenario, after 
peak the thin shell has been overtaken by the ejecta, and the 
luminosity is powered by the interaction of the ejecta with 
the slower, older wind. To explain this, the wind on day 711 
requires M w 10 -2 Mq yr -1 and v w ftt 100 km s _1 , which 
requires LBV-like mass loss (see above). The wind-blown 
bubble scenario therefore seems to unnecessarily complicate 
the scenario discussed above, which is needed to explain the 
late-time observations, as it requires the transition from a 
LBV to a Wolf-Rayet star. Yet the addition of this transition 
does not clearly model any observed properties that could 
not otherwise be explained. 

Additionally, we see no evidence for a normal SN in 
the years leading up to the initial detection of SN 2008iy 
as would be expected in the wind-blown bubble scenario. 
While it would be impossible to completely rule out a SN 
prior to the first detection of SN 2008iy on 2007 Sep. 13 (for 



instance, the field is virtually unobservable for ~3 months 
each year while it is behind the Sun), DS images dating to 

2004 show no evidence for such an event. These DS images 
have similar depths to those quoted in Section 2.1, which at 
the distance of the host galaxy corresponds to an absolute 
magnitude of Mi w —15.4 mag. This detection threshold is 
well below the mean peak magnitude for stripped-envelope 
SNe: Mr = -17.01, - 16.38, and -17 04 mag for SNe lb, Ic, 
and lib, respectively ( Li et al. |2010a |. Stronger constraints 
on the presence of a stripped-envelope SN prior to 2007 Sep. 
13 are based on the individual DS images. There are DS 
images of the field of SN 2008iy from 2004 Apr. and Jul., 

2005 May, Jun., and Jul., 2006 May and Jun., and 2007 Apr., 
May, Jun., and Jul. Adopting the mean peak magnitudes 
and mean light curves for stripped-envelope SNe from [0] 

(2010a I, we find that on average (note that there is 



et al. 



a great deal of dispersion in the peak absolute magnitude 
and duration of stripped-envelope SNe), SNe lb and lib are 
observable for ~70 days above the DS detection threshold, 
while SNe Ic are observable for only ~45 days above the 
DS detection threshold. In addition to not being detected 
during any of the months listed above, we can also rule out 
SNe lb and lib in the ~2 months prior to those listed above, 
while SNe Ic can be ruled out in the ~1.5 months prior to 
those above. Therefore, in the ~3.5 years prior to SN 2008iy 
we can rule out SNe lb or lib in ~21 out of 42 months, while 
SNe Ic can be ruled out for ~18 out of 42 months. 

Perhaps the most intriguing points about the progeni- 
tor of SN 2008iy are the implications for the progenitors of 
the very long-lived SNe Iln discussed in Section 4.2, such 
as SN 1988Z. It is interesting that episodic, LBV-like erup- 
tions, which we argue occurred shortly before SN 2008iy, 
may be applicable to SNe 1988Z, 1995N, and 1986J, each of 
which has relatively poor constraints on the actual date of 
explosion. 



SN 1988Z was first observed after peak (Stathakis & 



Sadler 19911, with an observational constraint on the rise 
time of < 250 d ays (|Turatto et al. |1993| . It was argued by 



Turatto et al. ( 1993 1 that the rise time of SN 1988Z must 



have been short, in part because at the time no other known 
SN looked quite like SN 1988Z. SN 2008iy, which is a virtual 
clone of SN 1988Z at late times, does provide an example of a 
SN with a long rise time and a peak absolute magnitude that 
is brighter than the discovery magnitude of SN 1988Z. This 
negates the necessity of a short rise to peak for SN 1988Z. We 
submit that it is possible that SN 1988Z had a long rise time, 
though we note that with no observations during the 250 
days prior to discovery, there is no way to definitely prove 
this claim one way or another. A long rise to peak would, 
however, further strengthen the similarity of SN 2008iy and 
SN 1988Z. 

Similarly, SNe 1986J and 1995N have poor constraints 
on the rise time. SN 1986J was discovered in the radio, with 
constraints on the explosion date ranging from 1982 to 1984; 
constraints on the optical rise are even worse, owing in part 
to the large extinction in the host galaxy, Ay ~ 2 mag (see 
Rupen et al. 1987 1. Typically, ~10 months prior to discovery 



is adopted as the explosion date for SN 1995N (see |Fox et al.| 
2000). However, this age of 10 months is based on the Ha 
profile of SN 1995N, which looked similar to the Ha profile 
of SN 1993N 10 months after SN 199 3N had been discovered 
( |Benetti, Bouchet fc Schwarz|[l995l ). We note that SNe Iln 
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are a very heterogeneous subclass, and in our experience 
spectral ages of SNe Iln are not reliable, especially at times 
greater than a few months. Thus, poor optical constraints 
mean that SNe 1986J, 1988Z, and 1995N all may have had 
long rise times, by the standards of typical SNe II whose rise 
times are < 1 week ( Li et al. |2010a |. 



4.4 Similarities to Very Luminous SNe Iln 

The long rise time and broad, symmetric evolution of 
SN 2008iy around peak optical emission is reminiscent of 



the Type Iln VLSNe 2006gy and 2008 fz ( |Ofek et al. [2007 
Smith et al. |2007| [Drake et al. |2009b|>. A pplying the "shell 
shock" model of Smith & McCray (20071, which was origi- 



nally developed to explain the light curve of SN 2006gy, does 
not provide a physically plausible scenario for SN 2008iy be- 
cause the relatively low M and large radius at peak imply 
that the CSM shell was not opaque. The "shell shock" model 
argues that the extreme luminosity of these SNe is powered 
via a shock running into a dense, optically thick shell. Ra- 
diation from the shock is thermalised by the optically thick 
gas and must diffuse out of the shell. 

Nevertheless, SN 2008iy does share a few similarities 
with the VLSNe Iln. Like SNe 2006gy and 2006tf, SN 2008iy 
shows evidence for a ~100 km s _1 pre-shock CSM wind (see 
Smith et aL~||2007[ |2008a[), and like SNe 2006gy, 2006tf, 



and 2008fz, SN 2008iy experienced a period of enhanced 
mass loss during the decades prior to explosion. The differ- 
ence between these systems is that the mass-loss rate from 
SN 2008iy, M « 0.01 M yr -1 , was less extreme than that 
from SNe 2006gy, 2006t f, 2008fz, where M « 1.0 M Q yr -1 
(e.g., Smith et al. |20f0 ). Also, the time between the period 
of enhanced mass loss and the SN was longer for SN 2008iy, 
~ 55 years, than fo r SNe 2006gy and 2006tf, < 10 years 
( Smith et al. | |2010 



2008a) 



The result of these differences 
is that SN 2008iy did not reach as extreme a peak luminos- 
ity and took considerably more time to reach peak optical 
output than the VLSNe Iln. Finally, if the late-time NIR 
excess is the result of an IR dust echo, this would be similar 



to SN 2006gy (Smith et al. ||2008b 


Agnoletto et al. ||2009 


Miller et al. 2009b I. Smith et al. 


( 2008b 1 and Miller et 


al. ( 2009b 1 argue that a giant eruption >1000 year prior 



to SN 2006gy could potentially create the dust shell giving 
rise to the echoes, and a similar shell around SN 2008iy may 
explain the observed NIR excess. 

The connection between VLSNe and very massive stars, 



such as r\ Car, was first made by Smith et al. ( 2007 1 because 
the large mass-loss rates needed to explain the extreme op- 
tical luminosity are reminiscent of r\ Car during the great 
eruption of 1843. It is interesting to note that were rj Car to 
explode as a SN today, its appearance would be less like the 
VLSNe Iln discussed above and more similar to SN 2008iy. 
The reason for this is that with the eruption happening >150 
yr ago, the ejected shell surrounding r\ Car has had sufficient 
time to expand (R « few x 10 17 cm) and become clumpy 
( Morse et al. |1998 1, with an optical depth that is now near 



unity (Davidson et al. 20011. Thus, r\ Car does not have 

cm), very optically thick shell 



SN. Instead, as the SN ejecta sweep through the large shell[^] 
they will successively overtake more and more mass within 
the clumpy CSM. This will likely give rise to a long rise time, 
similar to what was seen in SN 2008iy. 



4.5 The Host of SN 2008iy 

At the adopted distance to SN 2008iy, our detection of the 
host galaxy in stacked SDSS images (see Section 2.3) means 
that the host has an absolute magnitude of M r « —13.7 
mag. This is significantly fainter than the Small Magellanic 
Cloud (with Mv = —16.9 mag). We do not detect narrow 
emission lines from the host in our SN spectra in order to 
make a direct measurement of the metallicity. The prospects 
for such a measurement in the near future are not good: if 
SN 2008iy continues to evolve like SN 1988Z, light from the 
SN may dominate over light from the host galaxy for years 
to come. 

We can, however, estimate the metallicity based on 



the absolute magnitude of the host. Lee et al. (20061 use 
27 dwarf irregular galaxies to determine the luminosity- 
metallicity relation for low-mass galaxies, and following from 
their Equation 1, we find 12 + log(0/H) ~ 7.7, assuming 
Mb « M r . The uncertainty in this value is large, perhaps 
as great as ± 0.3 dex, based both on our assumption of 
B — r = mag for the host and the substantial scatter in 
the luminosity-metallicity relationship for faint galaxies ( |Lee| 
et al. |2006 find a scatter of 0.16 dex), yet this nevertheless 
shows that SN 2008iy occurred in a metal-poor galaxy. 

It has recently been suggested that the unusual SNe 
discovered by the non-targeted transient surveys, such as 
CRTS, preferentially occur in subluminous, possibly metal- 
poor host galaxies ( jMiller et al. |2009a| |Drake et al. |2009b I. 
SN 2008iy is yet another example of an unusual SN in a 
low-luminosity host. We caution, however, that a number of 
biases may be skewing initial impressions about the hosts 
of these unusual SNe. Preliminary calculations show that 



2009a Quimby et al. 



the compact {R < few x 10 10 
required to generate an extreme peak luminosity following a 



these VLSNe are rare (Miller et al. 

|| |2009[ ). Both the large peak luminosity and slow evolution 
of SN 2008iy mean that similar SNe would be easily de- 
tectable in the galaxies of targeted SN surveys. The lack of 
other examples of SNe with ~400 day rise times suggests 
that SN 2008iy is also rare. The CRTS, however, does not 
employ image subtraction during their search for transients; 
instead they use aperture photometry to find new sources, 
or sources with large increases in flux, and flag those as pos- 
sible transients. Consequently, their survey is biased toward 
the discovery of intrinsically bright SNe in faint host galax- 
ies. From the Lick Observatory Supernova Search (LOSS), 
we know that SNe Iln are rare, regardless of whether they 
are very luminous or have long rise times ( |Li et al. |2010b[ ). 
Yet, even LOSS, which targets specific galaxies, may be bi- 
ased in that it observes relatively few subluminous galaxies. 
New and upcoming surveys, such as the Palomar Transient 
Factory ( Law et al. ||2009 Rau et al~||2009 |, which employ 
image subtraction and survey large, non-targeted fields of 



8 The shell around r) Car is much bigger than the one around 
SN 2008 iy because th e wind speeds are greater, v WtV c a r = 600 
km s~ x (Smith 2006), and the nebula is older, ~165 yr relative 
to ~55 yr. 
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view, will not suffer from the same biases discussed above, 
and hence will be in a better position to address whether 
these unusual SNe preferentially occur in subluminous host 
galaxies. 



5 CONCLUSIONS 

We have reported on observations of the Type Iln SN 2008iy, 
which took ~400 days to reach peak optical output. There 
are few known SNe with optical rise times > 50 days, and 
SN 2008iy is the first with a rise time > 1 year. We argue 
that this long rise to peak is caused by the interaction of 
the SN ejecta with a dense CSM; radioactivity is unlikely 
to drive a 400-day rise. Furthermore, the late-time optical 
decay, ~0.4 mag (100 days) -1 , is slower than that of 56 Co, 
which provides further evidence that radioactive heating is 
not a dominant energy source for SN 2008iy. Spectroscopi- 
cally, SN 2008iy is very similar to SN 1988Z at late times. 
SN 1988Z is understood to have exploded in a dense, clumpy 
CSM (IChugai & Danziger||1994l), which we argue was also 



the case for SN 2008iy. We detect SN 2008iy in X-rays with 
a total luminosity Lx = (3.7 ± 1.2) x 10 41 erg s" 1 , which is 



similar to that of the Type Iln SNe 1988Z and 1995N (Fox 



et al. 2000). Similar to other SNe Iln, SN 2008iy has a grow- 



ing NIR-excess at late times. SN 2008iy had a peak absolute 
magnitude of M r « —19.1 mag and a total radiated energy 
of ~ 2 x 10 50 erg in the optical (assuming no bolometric 
correction) . 

The steady increase in optical luminosity over a ~400- 
day period means that the wind-density parameter, w, in- 
creased over a distance of ~1.7 xlO 16 cm from the SN. We 
propose two possible scenarios to explain this increase in 
w. (i) the progenitor experienced an episode of LBV-like, 
eruptive mass-loss ~55 years prior to the SN, or (ii) the 
wind speed of the progenitor was increasing during the years 
leading up to core collapse. We prefer the former scenario, as 
the later adds unnecessarily complicated pieces to the puzzle 
without providing a unique solution. Our favoured scenario 
provides yet another piece of evidence that some SNe Iln are 



connected to LBV-like progenitors (see Gal- Yam et al. 2007 



and references therein). We find that the host of SN 2008iy 
is a subluminous dwarf galaxy, though we caution against 
premature conclusions that unusual SNe, specifically those 
that are very luminous or have very long rise times, prefer- 
entially occur in low-mass dwarf galaxies. 

Finally, we close with some predictions for the late-time 
behaviour of SN 2008iy. There are examples of SNe Iln 
whose luminosity dramatically drops after the SN ejecta 
overtake the dense CSM (e.g., SN 1994W; |Sollerman" 



Cumming & Lundqvist 1998). However, the similarities to 



SN 1988Z suggest that SN 2008iy could continue interact- 
ing, and thus remain luminous, for several years. This would 
allow long-term monitoring in the radio, X-ray, and optical, 
like SN 1988Z (|Williams et al. |2002| jSchlegel fc Petre|2006| 



Aretxaga et al. ||1999[ ). We predict that SN 2008iy is a lu - 
minous radio source, like SN 1988Z ( |Van Dyk et al. |1993[ |, 
though we note that at z = 0.0411 deep observations may be 
necessary for detection. The increasing if s -band luminosity 
from day 560 to 710 is likely due to the presence of dust, 
and we predict that SN 2008iy will also be luminous in the 
mid-IR (~3— 5 /im). The data are insufficient to distinguish 



between newly formed dust or dust that was present prior to 
the SN, but future medium- to high-resolution spectroscopy 
could distinguish between these two cases. If new dust is be- 
ing formed, it should result in a systematic blueshift in the 
line profiles, as the dust creates an optically thick barrier to 
radiation from the receding SN ejecta. 
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